Many-Body Electronic Structure of Americium metal 
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We report computer based simulations of energetics, spectroscopy and electron-phonon interaction 
of americium using a novel spectral density functional method. This approach gives rise to a new 
concept of a many-body electronic structure and reveals the unexpected mixed valence regime of Am 
5f^ electrons which under pressure acquire the 5f^ valence state. This explains unique properties of 
Am and addresses the fundamental issue of how the localization delocalization edge is approached 
from the localized side in a closed shell system. 

PACS numbers: 



Artificially produced from Plutonium-239 in 1944, and 
widely used in smoke detectors Americium is the first 
transuranic actinide where 5f^ electrons become local- 
ized and form a closed relativistic subshell. Its recent 
high- pressure studies have drawn much attention as 
understanding volume behavior in actinides systems has 
important consequences on their storage and disposal. 
They have revealed that Am undergoes a series of struc- 
tural phase transitions (denoted hereafter as I, II, III, 
and IV) and reproduces at least two of the structures 
of another mysterious element, Plutonium, which links 
the physical behavior of all actinides materials to our 
fundamental understanding of bonding between their 5f- 
electrons. At ambient pressure Am I behaves as an or- 
dinary metal with slightly enhanced electrical resistivity 
/9(T=300K)=68 /ifi X cm and no sign of ordered or dis- 
ordered magnetism. This is standardly understood as a 
manifestation of ^F^ ground state singlet of 5f^ atomic 
configuration. However, the resistivity of Am raises al- 
most an order of magnitude and reaches its value of 500 
fiVl X cm at the orthorhombic structure of Am IV which is 
realized at pressures P above 16 GPa. The most promi- 
nent feature of the pressure P vs. volume V behavior is 
the existence of two distinct phases: the "soft" one which 
occurs in Am I through III as well as another "hard" 
phase realized in Am IV. On top of that a superconduc- 
tivity in Am was first predicted^ and then discovered^ 
with Tc raising from 0.5K in Am I to 2.2K in Am II, 
falling slightly in Am III and then exhibiting a sharp 
maximum in phase IV 0|. 

Understanding this unique behavior is a fundamental 
challenge in searching for a unified theory of actinides 
as the pressure driven delocalization of electrons is ap- 
proached here from the localized side which is very dif- 
ferent from Pu where originally delocalized electrons be- 
come localized with increasing volume. Thus, simple 
model Hamiltonians which contain qualitative features to 
produce complex energy landscapes with multiple solu- 
tions in open shell systems cannot be employed for stud- 
ies of closed shell materials: without incorporating re- 
alistic structures in the calculation, there is no hint of 



bistability in the model Hamiltonian approach. 

To address these issues in this work we introduce a 
novel many-body electronic structure method which al- 
lows us to uncover the physics of Am. It is based on 
dynamical mean field theory (DMFT), a modern many 
body technique for treating strongly correlated electronic 
systems in a non-perturbative manner |^ Q and at the 
same time has computational efficiency comparable with 
ordinary electronic structure calculations thus allowing 
us to deal with complicated crystal structures of real 
solids by self-consistent many-body calculations. Our 
new method considers the local Green function Gioc{'^) 
as a variable in the total energy functional and can be 
viewed as spectral density functional theory,^ ,^|^J. The 
advantage of such formulation as compared to original 
density functional theory^^ is a simultaneous access to 
energetics and local excitation spectra of materials with 
arbitrary strength of the local Coulomb interaction U. 

DMFT based spectral density functional approach re- 
quires self-consistent solutions of the Dyson equations 



[uj - iJo(k) - S](w)]G(k, w) = 1 



(1) 



for the one-electron Green function G'(k, w). The poles 
of its momentum integrated Gioc{^) contain information 
of the true local spectra of excitations. Here -ffo(k) is 
the effective one-electron Hamiltonian while S(ijj) is a lo- 
cal self-energy operator whose energy dependence makes 
the solution computationally very expensive. This so 
far has restricted applications of this promising many- 
body approach either to non-self-consistcnt determina- 
tions of spectraill| or to materials with simple crystal 
structuresliimTal. 

Our new approach greatly improves the speed of the 
calculation by recognizing that a signature of strong cor- 
relation effect results in appearance of several distinct 
features or satellites in the excitation spectrum. The ex- 
act self-energy of an interacting system can always be 
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represented by a pole expansion of the form 

I]M=E(oo)+5]-^ (2) 

I 

Remarkably, that such form of the self-energy allows us to 
I 



which relates our original matrix inversion required to 
find G(k, uj) (first element in the matrix from the right) 
to the matrix inversion in the extended "pole space" . 

The key insight is that the above form of the self energy 
with a few poles captures |^ all the central features 
of a correlated system and is an excellent approximation 
to the Greens function of the system for the purposes of 
obtaining the total energy. This has an important im- 
plication for the calculation of the electronic structure of 
the strongly correlated material: once pole expansion of 
the self-energy is established, the spectral density func- 
tional theory reduces to solving a "Kohn Sham-like" sys- 
tem of equations in an augmented space. The eigenstates 
here describe major atomic multiplet transitions as well 
as delocalized parts of the electronic states by separate 
auxiliary wave functions. Each wave function is not nor- 
malized to unity since it describes only part of the spec- 
tral weight for the electron leaving in the vicinity of a 
given energy, however the integral spectral weight over 
all energies is correctly normalized to one. While a self 
energy with a small number of poles may not capture 
subtle physics of damping important for incoherent ex- 
citations, the method gives us directly the dispersions 
of these spectral features, which are measurable in angle 
resolved photoemission. 

Thus, the concept of the electronic structure is general- 
ized to a strongly correlated situation. It ideally suits the 
description of such subtle regime as the proximity to the 
Mott transition where atomic multiplet structure appears 
simultaneously with strongly rcnormalizcd quasiparticle 
bands, a regime where traditional electronic structure 
methods fail. 

Here, we study the properties of Am under pressure us- 
ing this newly implemented matrix expansion algorithm 
for spectral density functional calculations within a full 
potential version of the linear mufhn-tin orbital method 
|1?tI |. In this approach, the s,p,d electrons are assumed to 
be weakly correlated and well described within such pop- 
ular approximations to the density functional theory as 
the local density approximation (LDA) including gradi- 



replace the non-linear (over energy) Dyson equation by 
a linear Schroedinger-like equation in extended subset 
of "pole states". This is clear due to a mathematical 
identity 



(3) 

I 

ent corrections (GGA). The correlated f electrons require 
dynamical treatment using DMFT. Both, the spin-orbit 
as well as the Hund's couplings are competing in Ameri- 
cium and need to be taken into account. The former is a 
one-body term and enters through the LDA Hamiltonian, 
while the second is contained in the local Coulomb repul- 
sion, which is conveniently expressed via Slater constants 
F'^*^ The value of the most important term f(°)=U is 
around 4.5 eV which is suggested from various atomic 
spectroscopy data and our previous studies of Plutonium. 
For the remaining constants we take the atomic values 
F(2)=7.2 eV, FW=4.8 eV, F(6)=3.6 eY^. 

For the purpose of the total-energy calculation, the f- 
electron self-energy is approximated by its atomic value 
which is obtained by the exact diagonalization technique. 
This is known as the Hubbard I approximation. The 
probabilities to find the f-shcU in its given many body 
state are directly accessed within this method and give 
us the insight into valence of the material. 

Our calculation reproduces the well known fact that 
the f electrons in Am at zero pressure exists in a f^ yF^ 
configuration. This is illustrated in Fig. 1 by plotting 
density of states and energy bands reflecting the atomic 
multiplet transitions which demonstrate our novel ma- 
trix expansion algorithm. Our calculated one-electron 
density of states shown on the right consists of several 
distinct features related to ^f^ electron removal and 

— s-f'' electron addition processes. One can see that oc- 
cupied part of the spectrum is well compared with the 
available photoemission experimentsllTl . thus advancing 
previous bulk-surface interpretation|lS| as well as den- 
sity functional based calculations [iSl |2(J. The idea of 
our new method is to model three major satellites re- 
lated to 7F° ^ eH^/', 7F° ^ 8S^/2 and yF^ ^ qP^^ 
transitions with two pole self-energy and resolve them 
as many-body energy bands. This is illustrated on left 
part of Fig. 1 where the method is seen to capture the 
spectral weight related to electron removal and addition 
processes by a set of eigenstates located near -3 eV bind- 
ing energy (blue) and by two sets of eigenstates located 
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FIG. 1: Calculated many body set of energy bands (left) 
which models the one-electron density of states (right) of 
Americium metal. Experimental photoemission spectrumfr3| 
is also shown by open circles. 



at +1 (green) and +3 eV (red). A simplified fee struc- 
ture with equilibrium atomie volume of Am I was used 
to generate the data in Fig 1. 

The computational speed gained by this algorithm al- 
lows us to study complicated crystal structures of Am. 
In particular, the existence of soft and hard phases in its 
equation of state can be predicted via our self-consistent 
total energy calculations. This is illustrated in Fig. 2 
where P(V) behavior reconstructed from the total en- 
ergy data of phases I through IV is plotted and com- 
pared with the recent experiment. For Am I we predict 
the equilibrium volume equal to 27.4 A'^/atom which is 
only 7% less than the experiment together with the bulk 
modulus equal to 450 kBar close to the experimentally 
deduced values lying within 400-450 kBar. The pressure 
ranges of all other structures are correctly reproduced. 
A compressibility of highly pressurized Am IV structure 
is found to behave similarly to experimentally observed 
"hard" phase which indicates that f-electrons start par- 
ticipating in bonding. Some discrepancy between the cal- 
culated and the measured compressibility of Am III can 
be observed in Fig. 2. This is likely to be due to sim- 
plified impurity solver or due to the uncertainties in the 
estimates of the Hubbard U. 

To gain theoretical insight and understand the origin of 
localization-delocalization transition we now discuss the 
behavior of the electronic structure under pressure. To 
see how the increase of hybridization among f-electrons 
affects the physical properties of Am, we carry out sub- 
sequent refined calculations, by replacing the Hubbard I 
approximation by a more precise one-crossing approxi- 
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FIG. 2: Calculated equation of state for various crystal struc- 
tures of Am metal. Experimental datal are shown by solid 
cicrles. 



mation (OCA) method|13ll22| to solve the Anderson im- 
purity problem. Due to numerical complexity of the ap- 
proach, we calculate only spectral functions using the fee 
structure of Am and omit self-consistent determination 
of the energy. 

Fig. 3 shows the density of states for Am at three 
different volumes V=Vo,V=0.76Vo and V=0.63Vo which 
covers the lattice spacing of the entire phase diagram dis- 
cussed above. Upon compression, the remarkable effect is 
observed as peak near the Fermi level gets pushed down 
while a resonance (small shoulder) starts forming at E/ 
and becomes more pronounced with increasing pressure. 
The f^ ground state of the atom starts admixing an f'' 
configuration with a very large total spin of J=7/2. Due 
to hybridization with the spd bands, this large spin gets 
screened thus lowering the energy of the system. This is 
the famous Kondo mechanism, and the energy gain in- 
creases as the hybridization increases by applying pres- 
sure. 

The admixture of the f configuration is counterintu- 
itive. Naively one expects that application of pressure 
results in lowering the Fermi level in the spd band (which 
contains only 3 electrons) which then moves towards the 
f level. This reduces the occupancy of the f, admixing an 
f^ configuration, an effect that is known to induce mixed 
valence in Sm compounds |2^ . Our first principles cal- 
culations reveal that while the position of the bare f level 
in Am indeed moves upwards relative to the Fermi level, 
the energy to absorb an electron to reach f configuration 
is much smaller than the energy to remove an f electron 
and transfer it to the Fermi level hence reaching the f^ 
configuration. Application of pressure reduces the en- 
ergy of the f^ configuration in the bath of spd electrons 
by the gain in hybridization, and the resulting energy 
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FIG. 3: Pressure dependence of Am density of states calcu- 
lated using the one crossing approximation method. 



gain is sufficient to compensate the increase in the dis- 
tance between the bare f level and the Fermi level. In 
this way, the valence is also changed from 3 to approx- 
imately 2.8 under extreme pressure. This conffi'ms the 
attribution of the rise in resistivity to mixed valence and 
provides the microscopic mechanism of this phenomenon. 
At normal pressure, the f shell is essentially closed and 
unable to scatter the conducting bands. This results in 
small resistivity. Under pressure, the reaction i* + spd 
—> f^ becomes more energetically allowed resulting in the 
growth of resistivity. 

It is interesting to compare the DMFT predictions for 
Am with the predictions of the density functional theory. 
The non-magnetic GGA calculation falls catastrophically 
in reproducing the theoretical equilibrium volume of the 
soft phase and underestimates it by about 50%. When 
spin polarization is allowed, the GGA eventually recovers 
most of this error but converges to the wrong magnetic 
state with its total (spin p lus orbital) moment of about 
6 Bohr magneton |l9t |2(J. This prediction is at odds 
with the experimentally established f^ ground state sin- 
glet rF^ ■ Similar findings have been reported when ap- 
plying a disordered local moment method [2l| . This error 
is the result of neglecting an important correlation eff'ect. 
The Kohn Sham spectrum of Am describes an f level with 
a small spin orbit splitting between ^^'^ and f^'^ (of the 
order of 1 eV) which leads to two energy bands located 
just near and which are unstable against magnetism. 
In Am the Coulomb interaction increases this splitting 
by the value of U (of the order of 4 eV) and leads to 
the atomic-like yF^ many-body state. This stabilization 
of the closed shell due to the Hubbard interaction U is 
absent in the band calculation. 

We finally estimate the superconducting critical tem- 



perature by computing from first principles [2J| the 
electron-phonon coupling of the electrons in the presence 
of correlations. For this purpose we have extended a 
newly developed dynamical mean field based linear re- 
sponse method, which has previously proven to provide 
accurate phonon spectra in correlated systems . 
We estimate the coupling constant which comes out to 
be sufficiently high (^0.5) to predict superconductivity 
of the order of 1 K. The occurrence of the first maximum 
in experimental Tc vs pressure dependence, can then be 
understood as the result of the variation of the spd den- 
sity of states which first increases as a result of a band 
structure effect but then eventually decreases as the hy- 
bridization with the f electron grows with the increase of 
mixed valence. 

To summarize, here we provided a first order picture 
of electronic properties of Am metal but the necessity 
of its further studies is apparent. These require exten- 
sions of our methods to evaluate the electron-phonon and 
the Coulomb interactions among quasiparticles in a full 
fledged mixed valence state and will be carried out in the 
future work. 
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